Strongly Correlated States of Ultracold Rotating Dipolar Fermi Gases 
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We study strongly correlated ground and excited states of rotating quasi-2D Fermi gases constituted of a 
small number of dipole-dipole interacting particles with dipole moments polarized perpendicular to the plane 
of motion. As the number of atoms grows, the system enters an intermediate regime, where ground states are 
subject to a competition between distinct bulk-edge configurations. This effect obscures their description in 
terms of composite fermions and leads to the appearance of novel composite fermion quasi-hole states. In the 
presence of dipolar interactions, the principal Laughlin state at filling v — 1/3 exhibits a substantial energy gap 
for neutral (total angular momentum conserving) excitations, and is well-described as an incompressible Fermi 
liquid. Instead, at lower fillings, the ground state structure favors crystalline order. 
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Some of the most fascinating challenges of modern atomic 
and molecular physics arguably concern ultracold dipolar 
quantum gases |Q]]. The recent experimental realization of a 
quantum degenerate dipolar Bose gas of Chromium [2], and 
the progress in trapping and cooling of dipolar molecules ^ 
have opened the path towards ultracold quantum gases with 
dominant dipole interactions. Particularly interesting in this 
context are rotating dipolar gases (RDG). Bose-Einstein con- 
densates of RDGs exhibit novel forms of vortex lattices, e. g., 
square, stripe- and bubble- "crystal" lattices |4|]. The stability 
of these phases in the lowest Landau level was recently in- 
vestigated L5]. We have demonstrated that the quasi-hole gap 
survives the large N limit for fermionic RDGs [6] . This prop- 
erty makes them perfect candidates to approach the strongly 
correlated regime, and to realize Laughlin liquids (cf. |7|]) at 
filling v = 1/3, and quantum Wigner crystals at v < 1/7 Js|] 
for a mesoscopic number of atoms N ~ 50 — 100. Lately, 
Rezayi et al. L9Q have shown that the presence of a small 
amount of dipole-dipole interactions stabilizes the so-called 
bosonic Rezayi-Read state at v = 3/2 whose excitations are 
both fractional and non-Abelian. 

In this letter, we investigate quasi-2D microscopic and 
mesoscopic clouds of fermions whose dipole moments are po- 
larized perpendicular to their plane of motion. Results are 
obtained by exact diagonalization of the Hamiltonian of the 
system for N = 3 — 12 particles. We study the lowest Lan- 
dau level (LLL) physics, in particular the crossover from the 
weakly interacting Fermi liquid regime to strongly correlated 
fractional quantum Hall-like states. As N grows, the ground 
state particles either organize themselves in the bulk or crys- 
tallize on the edge. This competition is based on inter-particle 
correlations and exceeds the scope of the effective theory of 
composite fermions ifioll . It leads to the appearance of novel 
ground states with a quasi-hole at the origin which are linked 
to unfavored composite-fermion states. The analysis of the 
principal Laughlin state at filling v — 1/3 reveals that its low 
energy excitations in the angular momentum subspaces where 
quasi-holes are to be found, correspond to edge excitations 



rather than topological quasi-hole defects. Instead, "neutral" 
particle-hole pair excitations are already substantially gapped 
in these microscopic systems. 

We consider a system of TV dipolar fermions rotating in an 
axially symmetric harmonic trapping potential strongly con- 
fined in the direction of the axis of rotation. Along this z- 
axis, the dipole moments, as well as spins are assumed to be 
aligned. Various ways of experimental realizations of ultra- 
cold dipolar gases are discussed in HI]. In case of low tem- 
perature T and weak chemical potential \i with respect to the 
axial confinement ou z , the gas is effectively two-dimensional, 
and the Hamiltonian of the system in the rotating reference 
frame reads 
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Here, uq <C uj z is the radial trap frequency, Q is the fre- 
quency of rotation, M is the mass of the particles, Vd = 
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is the dipolar interaction potential (rotationally 



invariant with respect to the z-axis), d is the dipole moment, 
and Yj = xje x + yje y is the position vector of the j-th parti- 
cle. The first term of (Q} is formally equivalent to the Landau 
Hamiltonian of particles with mass M and charge e moving 
in a constant magnetic field of strength B = 2MVtc/e per- 
pendicular to their plane of motion. The spectrum of 7YLandau 
consists of equidistantly spaced and highly degenerate Lan- 
dau levels (LL) with energies e n = hw c {n + 1/2) where 
oo c = 20. We denote by TVll = 1/2tt1 2 the number of states 
per unit area in each LL, where I = ^h/muu c is the magnetic 
length. Given a homogeneous fermionic surface density rif, 
the filling factor v = 2irl 2 nf can be defined and denotes the 
fraction of occupied Landau levels. Even though the above 
definition applies to infinite homogeneous systems, it may be 
used for finite systems as a suitable truncation of the Hilbert 
space at specific angular momenta. The second term in ([T]) ac- 
counts for a rotationally induced effective reduction of the trap 
strength. For critical rotation O — > uo, the confining potential 
vanishes. 
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In the following, it is assumed that particles solely occupy 
the LLL lUlll . Restricted to the LLL and rewritten in second 
quantized form, the Hamiltonian (Q]) reads 

H = hw Q N + H(UJ0- £ z + Vl2M a mi a L 2 a m4am 3 (2) 

mi ,...,7714 

where N and are the total number and z-component an- 
gular momentum operators, . creates a particle with angu- 
lar momentum m^, and V1234 = \{™>i m 2 1 Vd 1^3 777,4) is the 
matrix element of interaction expressed in the Fock-Darwin 
single particle angular momentum basis II 1211 . Due to circu- 
lar symmetry, © can be diagonalized blockwise for a given 
L z . Calculations have been performed for N = 3 to 12 parti- 
cles with complete and Davidson block diagonalization tech- 
niques W13H . respectively. In the LLL Hamiltonian, besides 
the first constant term, there are two competing contributions 
to the energy: the kinetic term linear in L and the strength 
of interaction, which scales as d 2 /l 3 = 2huo(ad/l) with 
ad = Md 2 /h 2 . The natural unit of energy is Hujo, whereas 
distances are measured in I; from now on, we set ad/ 1 = 1. 
We analyze the ground state interaction energy for a given L z . 
It reveals plateaus, clearly visible for small 7V(see Fig.Q]). 

FIG. 1: Interaction contribution to the GS energy as a function of L z 
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Ground state candidates are the first states of these plateaus 
where a downward cusp occurs in the spectrum. By tuning 
the rotational frequency, some of these states are selected as 
true ground states at specific "magic" angular momenta as de- 
picted in Fig. [2l For relatively low values of Q, the ground 
state is the filled LL state at v = 1, which is completely in- 
sensitive to type and strength of interaction as long as the LLL 
approximation holds. 

When Q is continuously increased, the system evolves from 
the weakly interacting regime (where the second term of 
Eq.© is larger than the third) to strongly correlated states 
(where the interaction is dominant). In case of short-range in- 
teractions, this process terminates at Lg = v ~ 1 N(N — l)/2, 
where the fermionic (bosonic) Laughlin states at filling vo = 
1/3(1/2) become the true ground states. 

The existence of a final L z is due to the fact that the 
ground state contact interaction energy vanishes for L z > Lg. 



FIG. 2: Ground state angular momentum series over a = ouo — H 
(the divergence at a = is not shown). 
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The long range nature of dipolar interactions lifts this de- 
generacy. Thus, the whole principal series of fillings, i.e., 
v = l/(2m + 1) for fermionic gases, is accessible. To 
reveal the internal structures of relevant states, we consider 
the density-density correlation function p(r, fo), which rep- 
resents the conditional probability to find one atom at r when 
another is simultaneously located at fo , 

N 

p(r, f ) = ^2 ^( r i r)8(f* k - fo) . (3) 

j<k 

It is crucial to analyze second order correlations, since the GS 
density only reveals radially symmetric contributions. Fig. [3] 
depicts p(r, fo) for a selection of ground states with N — 10 
particles. Generally, for small N, starting from the unstruc- 

FIG. 3: Ground state density-density correlation functions p(r, fo) 
for N = 10 dipolar fermions at L z =(top) 45,80, (centered) 90,93, 
(bottom) 103, 117 with fo set to the maximum of the density, which 
occurs at the edge. 
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tured maximum-density-droplet at L z = N(N — l)/2 (top 
left) with v — 1, a fraction of atoms starts to arrange itself in 
a "crystal" on the edge leaving a residual "Fermi sea" at the 
center(top right), until the correlations are homogeneously es- 
tablished at Lg = 3N(N - l)/2. For larger N, however, 
the energetically favorable number of edge atoms in the crys- 
tal changes "irregularly" with L z (right column); the system 
seems to be "frustrated" in this respect. In this regime, the 
bulk has to reorganize itself accordingly, and novel fascinat- 
ing states with a hole at the center of the bulk appear (centered 
and bottom left). 

In order to understand the "magic" L z numbers, Jain et al 
iTToh have proposed to model the interacting system in terms 
of an effective theory of non-interacting composite fermions 
(CF) for electrons in a quantum dot. In this ansatz, each 
fermion captures an even number of quantum fluxes, and the 
wave function reads 

*CF({*i» = #AlL J H(Zj - ^) 2m ^Landau > • (4) 
[j<k ) 

Here, m denotes the number of flux pairs, ^Landau is an TV- 
particle eigenfunction of 7YLandau, and Vlll is the LLL pro- 



jector. Motivated by the success of CFs for fermions, Cooper 
and Wilkin have nicely adapted this scheme to bosonic contact 
interacting gases full . Following this idea, we compared the 
series of true GS for non-interacting composites and dipolar 
fermions. As long as there is no real bulk in the system, i. e., 
for N < 7, the ground state series nearly identically match 
with the predictions of the effective CF theory. Furthermore, 
overlap calculations for electrons and short-range interacting 
bosons ESQ suggest that similar results will hold for dipo- 
lar fermions. For bigger N, however, important deviations 
from CF theory, in particular in the intermediate regime, oc- 
cur (see, for instance, the CF and true "magic" numbers in 
the table [5]). This deficiency of the CF theory, already com- 
mented in Ref. 11411 . is more clearly related to the "frustra- 
tion" effect and reorganization of the bulk in the present case. 
For N = 10, ground states with a density defect at the cen- 
ter are found at L z = 90, 95 and 103. They either have no 
(L z = 95), or are hardly similar to the corresponding "magic" 
CF states. Instead, the analysis for the accessible range of sys- 
tems strongly suggests that these states originate from "par- 
ent" states which are boosted by N quanta of angular momen- 
tum, e. g., L z = 80, 85 and 93 for the above series. None of 
these parent states has its "magic" analogues in the CF model. 
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To a good approximation, the close connection between 
parent and boosted state (e.g. L = 93 and 103 respectively) 
can be understood as a quasi-hole excitation, which is analyt- 
ically represented by 

*qh(te}) = A/" | J[( Zj - Co) m ^ | , (6) 

where ("o = 0, and m = 1 for = 1/3. Indeed, the above 
wave function proves to be a good approximation for the states 
with L z = L z arent + N: even though the exact states do 
not reveal a true topological defect at the origin due to finite 
size, the density at the origin scales as 1/N and decreases 
from 0.14 to c± 0.10 as N varies from 7 to 11, while overlaps 
with the exact ground states grow from 0.6 to 0.7 despite the 
simultaneous significant increase of the relevant Hilbert space 
dimension. 

We have also studied in detail low energy excitations of the 
dipolar Laughlin state at v = 1/3, which for A 7 " > 10 con- 
sists of a significant bulk, surrounded by a practically "melted 
crystal" at the edge (see Fig. [4]). In order to relate our results 
to the ones obtained in the thermodynamical limit in Ref. 
it is necessary to identify the finite-size quasi-hole excitations 
in the spectrum. To map out irrelevant candidates, it is reason- 
able to appropriately truncate the Hilbert space by a finite size 



FIG. 4: Density-density correlation function p(r, fo) of the Laughlin 
state for A = 12 dipolar fermions with fo chosen at the maximum 
of the density, which occurs at the edge. 
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filling factor. We use the approach of Ref. II X 511 . and fix v by 
imposing the constraint m v < (TV — l)/v on the maximum 
single particle angular momentum in the TV-particle Fock ba- 
sis. If a quasi-hole (quasi-particle) is nucleated, the filling fac- 
tor is lowered (raised) by a correction of the order of 1/N, i. e., 
m v — > m v ± 1. This significantly reduces the basis dimension 
of the system and maps out irrelevant states. If the rotational 
frequency Q is tuned to favor ^y 3 as the true ground state, 
the lower lying excitations for Ly 3 < L z < Ly 3 + N turn 
out to be gapped by Ae 1/N. For short-range interacting 
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bosonic systems, these states have been partially identified as 
edge excitations that carry quanta of angular momentum Il6ll . 
To identify the excited states with a quasi-hole at the centre, 
the ground states at L z qh = L^ 3 + N are calculated with 
adapted m v + 1. These states closely follow the lowest en- 
ergy branch of edge excitations, and should be regarded as 
such, rather than quasi-holes. This is strongly supported by 
the fact that the substantial density dip for N = 3 continu- 
ously vanishes for increasing number of particles. Additional 
information about these edge excitations can be obtained by 
direct "naked-eye" inspection of Fig. HI For the "standard" 
Laughlin state, the edge excitations are described by a hydro- 
dynamical Luttinger liquid theory [17], and their amplitude de- 
cays as sin -3 ((/>), where <\> is the "angular distance" from ro 
along the edge; here the decay is much slower, which is yet 
another signature of the finite size effects. 

The above considerations imply that the only reliable quan- 
tity that remains to estimate the quasi-hole energy in the large 
N limit is the neutral gap at fixed L z . This gap remains more 
or less constant for smaller TVs, and graduately increases as 
the number of bulk particles in the system and TV grow. For 
the data of Ref. [6] (M = 30 a.m.u., d = 0.5 Debye, and a trap 
frequency of 2tt x 10 3 Hz), the gap is of the order of percents 
of 2Hujo, i. e., it is substantial, but an order of magnitude less 
compared to the value estimated in ((J in the large N limit. 
Obviously, the lack of true bulk behavior in the investigated 
mesoscopic samples is responsible for this effect. 

In the semiconductor fractional quantum Hall effect, 
Wigner crystals, i. e., specific charge-density waves, were 
discussed as competing ground states to Laughlin liquids. 
The interplay of quantum fluctuations with the interaction en- 
ergy have proven crystalline Wigner order to be favorable for 
low enough fillings for electrons. This behavior intuitively 
changes in the case of dipolar particles as the interaction en- 
ergy scales differently in the density of the particles. A very 
recent detailed analysis of this issue has confirmed the sta- 
bility estimates for Laughlin and Wigner states of Ref. (6|] 
for systems constituted of 50-200 particles |8[]. Fehrmann 
et al. have proven that phonon excitations destabilize the 
Wigner crystal for fillings v > 1/7, which melts into a Laugh- 
lin liquid. In the microsystems discussed in this letter, this 
crossover is strongly supported. For fillings v < 1/5, den- 
sity profiles significantly deviate from the Laughlin state and 
show a clear pattern of hexagonal order for N = 6 particles 
as depicted in Fig. [5] where the outer edge ring of the state is 
constituted of five particles. Even though speaking of a true 
"crystal" is certainly not applicable in such a system, signa- 
tures of crystalline order are clearly present. 

Summarizing, we have studied in detail ground and excited 
states of quasi-2D ultracold rotating dipolar Fermi gases. By 
exact diagonalization methods, we studied systems up to 12 
particles. We have identified novel kinds of strongly corre- 
lated states in the intermediate regime, i. e., pseudo-hole exci- 
tations of CF states. Calculation of the substantial gap in the 
excitation spectrum of the dipolar Laughlin state at v = 1/3 
prove the accessibility of fractional quantum Hall states in 



FIG. 5: Radial density p(r) of Laughlin states (circles) and true dipo- 
lar ground states (triangles) for N — 6 at different filling factors. 
Apart from five particles which constitute the edge ring, localization 
of the sixth particle at the origin is clearly visible for v < 1/5. 
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these microsystems. At lower fillings, interactions favor crys- 
talline order. Rotating dipolar gases are thus very suitable 
candidates to realize Laughlin-like and more exotic quantum 
liquids, as well as their crossover behavior to Wigner crystals. 
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